Coherent Dynamics of a Josephson Charge Qubit 
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We have fabricated a Josephson charge qubit by capacitively coupling a single-Cooper-pair box 
(SCB) to an electrometer based upon a single-electron transistor configured for radio-frequency 
readout (RF-SET). Charge quantization of 2e is observed and microwave spectroscopy is used to 
extract the Josephson and charging energies of the box. We perform coherent manipulation of the 
SCB by using very fast DC pulses and observe quantum oscillations in time of the charge that persist 
to ~ 10ns. The observed contrast of the oscillations is high and agrees with that expected from the 
finite Ej/Ec ratio and finite rise-time of the DC pulses. In addition, we are able to demonstrate 
nearly 100% initial charge state polarization. We also present a method to determine the relaxation 
time Ti when it is shorter than the measurement time T meaa . 
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Although a large number of physical systems have been 
suggested as potential implementations of qubits, solid 
state systems are attractive in that they offer a realis- 
tic possibility of scaling to a large number of interact- 
ing qubits. Recently there has been considerable exper- 
imental progress using superconducting microelectronic 
circuits to construct artificial two-level systems. A va- 
riety of relative Josephson and Coloumb energy scales 
have been used to construct qubits based upon a single- 
Cooper-pair box EI 0] and flux qubits based upon a 3- 
junction loop Coherence times of the order of 0.5 

microseconds have been achieved for a single-Cooper-pair 
box qubit 0]. Rabi oscillations between energy levels of 
a single large tunnel junction have also been observed 
[ELE1< Despite the encouraging results, one aspect that 
is not well understood concerns the contrast of the oscil- 
lations, which in all previously reported experiments is 
smaller than expected. 

The experimental systems reported so far can also be 
distinguished by the readout method and the manner 
of performing single qubit rotations. The first demon- 
stration of coherent control of a single-Cooper-pair box 
(SCB) |lj employed a weakly coupled probe junction to 
determine the charge on the island. In the more recent 
experiment reported by Vion et al. |2| , the SCB was incor- 
porated into a loop containing a large tunnel junction, for 
which the switching current depends on the state of the 
SCB. Switching current measurements of SQUIDS have 
also been used for flux and phase- type qubits 0, EL EL El < 
Nakamura et a/.Q performed single qubit rotations by 
applying very fast DC pulses to a gate lead in order to 
quickly move the SCB into, and away from the charge 
degeneracy point. This technique produces qubit rota- 
tions with an operation time that can be of the order 
of the natural oscillation period. Other experiments uti- 
lize microwave pulses to perform NMR-like rotations of 
the qubit @,0,ME|- The latter approach requires less 
stringent microwave engineering, since RF-rotations can 
be accomplished with pulses that are more than order 
of magnitude slower in rise time and duration than the 
natural oscillation period. By using fast pulses, however, 
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FIG. 1: (a) Scanning electron micrograph of a sample. The 
device consists of a Cooper-pair box and SET electrometer 
and was fabricated from an aluminum (lighter regions) evap- 
orated onto an oxidized Si substrate (darker regions), (b) 
Circuit diagram of box and electrometer. 



single qubit rotations can be performed approximately 
20 times faster than with RF-rotations. 

In this letter, we report measurements made on a SCB- 
type qubit with very fast DC pulses used to effect the 
qubit rotations, as in Nakamura et o/.0j. For our qubit, 
however, the readout system consists of a single-electron- 
transistor (SET) capacitively coupled to the SCB and 
configured for radio- frequency readout (RF-SET) . By in- 
corporating a SET into a tank circuit, RF-SET electrom- 
eters can be made that are both fast0 and sensitive H, 
and are well-suited for SCB-qubit readout^, ^J. One 
advantage of using RF-SET readout is that it can eas- 
ily be turned on and off, although for the measurements 
reported here it is operated in a continuous [Tof. More 
significanly, the RF-SET is fundamentally different from 
qubit readout based on swithing currents in that it in- 
volves a weak measurement of the charge in contrast to a 
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yes or no answer. This could be a significant advantage 
for readout of multiple qubit systems. For example, a 
switching current measurement on one qubit would also 
directly affect the other qubits due to its strong interac- 
tion. Furthermore, switching current measurements are 
inherently stochastic in time — an important considera- 
tion if one is interested in measuring time correlations. 

Electron beam lithography and double-angle shadow 
evaporation of aluminum films onto an oxidized silicon 
substrate were used to fabricate the combined SCB-SET 
system (see Fig. 1)0. The SCB box consists of a low- 
capacitance superconducting island connected to a super- 
conducting reservoir by two parallel junctions that define 
a low-inductance SQUID loop. An additional gate lead 
placed near the island is used to change the electrostatic 
potential of the island by a gate voltage V g through a 
gate capacitance C g . By adjusting the magnetic flux $ 
through the loop, the effective Josephson energy is tun- 
able as Ej = Sj la:E |cos(7r$/$o)| J where $ is the mag- 
netic flux quantum (h/2e). Coupling of the SCB to the 
SET was accomplished by extending part of the SET is- 
land to the proximity of the box island and resulted in 
a weak dimensionless coupling, Cj/Ce of 0.5 — 4% for 
samples with slightly different geometries and total SCB 
capacitances Cj> Cj is the capacitance between the SCB 
and SET islands. 

The samples were placed at the mixing chamber of 
a dilution refrigerator with a base temperature of ap- 
proximately 20mK. An external superconducting magnet 
was used to produce a large field parallel to the plane of 
the samples and a small superconducting coil close to 
the sample was used to produce a field perpendicular to 
the SQUID loop. This allowed us to independently sup- 
press the superconducting gap of the aluminum film, and 
change the effective Josephson energy of the Cooper-pair 
box. All control lines were filtered by a combination of 
low-pass and stainless steel and copper-powder filters. 
To present sharp pulses to the SCB gate, a high fre- 
quency coaxial line was used having lOdB attenuation 
at 4.2K and 20dB at 1.5K. The tank circuit had a res- 
onant frequency of 380MHz and the RF-SET electrom- 
eter was biased near a feature in the IV versus V g ET 
landscape known as the double Josephson quasi-particle 
peak (DJQP) [l2j, such that the bias current through the 
SET was typically 200-300pA. Under these conditions, 
the electrometer sensitivity for the combined SCB-SET 
system was found to be 50 xlO -6 e/\/Hz with a band- 
width of 10MHz. 

To make an artificial two-level system, the energy 
scales of the SCB must be chosen so that fcsT < 
Ec,Ej > A, where Ec is the energy cost required to 
add a single electron to the island and is set by total 
capacitance of the island Eq = e 2 /2Cs- If A is large 
enough compared to Ec, then the ground states will be 
even parity states that differ only by the average number 
of Cooper-pairs on the island. The effective Hamiltonian 
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FIG. 2: (a) Ground state and first excited state energies ver- 
sus n g for Ej /Ec = .7 (solid line) and Ej = (dotted line), 
(b) Coulomb staircases, Qbox versus n g for a sample with 
E c /h = 9.25 ± 0.10GHz and EJ ax /h = 20.24 ± 0.10GHz. In 
(i) a staircase taken from a sample under microwave radia- 
tion with a frequency of 35GHz is shown. In (ii), we show 
data where Ej/h has been suppressed to 11.4GHz by apply- 
ing a small magnetic field through the SQUID loop. The 
dotted lines in (i) and (ii) show the calculated Qbox vs n g for 
the respective Josephson energies, (iii) Measured Qbox vs n g 
(solid line) when a fast pulse train is applied to the SCB gate 
having an amplitude 0.88e, width At = 130ps and repetition 
time T R = 130ns (here Ej/h = 4GHz). The dashed line 
is the measured Qbox without the pulse train, (c) Avoided 
level crossing shown by the positions of the spectral peaks 
and dips, An g — (n g hp — n g eak )/2, for varying microwave fre- 
quency, (d) Josephson energy Ej versus applied flux $/$o 
through the SQUID loop determined from spectroscopic data 
(circles) and coherent oscillation data (crosses). This sample 
had E c /h = 17.45 ± 0.13GHz and EJ ax /h = 15.5 ± 0.2GHz. 

of the box, including the Josephson coupling, is given by 
// = 4S c ^(n- % ) 2 |n)(n|-^^(|n+l)(n| + |n)(n+l|), 

n n 

where we define gate charge as n g — C g V g /2e—no, and uq 
is the offset charge due to stray charges near the box. Fig. 
2(a) shows the ground state and first excited state energy 
bands calculated using the Hamiltonian in Eqn.Q. On 
account of the 2e periodicity of the system, we limit the 
discussion to gate charge < n g < 1. 

In thermodynamic equilibrium, the actual quantity 
that determines the island parity is the even-odd free 
energy difference A(T), which differs from A due to 
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entropic considerations. 13J. Measurements were made 
for several different samples with Coulomb enegies 
E c /k B {E c /h) ranging from 0.43 to l.QSK (9-34GHz) 
and Ej/Eq ratios of 0.4 to 2.25. Coloumb staircases were 
measured by slowly ramping the box gate charge (at a 
rate of ~ 40e/s), and measuring the power reflected from 
the RF-SET tank circuit. Examples of the measured box 
charge Qbox{n g ) is shown in Fig. 2b. Our samples show 
2e periodicity when E c /ftB < IK, well below the gap 
A/ks = 2.2K. For samples with larger E c , the staircases 
acquire an extra step around n g — 0.5 due to poisoning 
by non-equilibrium quasiparticles. 

The characteristic energies Ec and E r j ax of the 
SCB can be determined directly using microwave 
spectroscopy When monochromatic microwaves are 
applied to the SCB gate, resonant peaks and dips oc- 
cur in the measured Qbox vs n g staircase when the mi- 
crowave photon energy matches the energy level splitting 
Voi (see Fig. 2(b)). As the microwave frequency is re- 
duced, the positions of these peaks exhibit an avoided 
level crossing. Within the two-level approximation, the 
positions of the single-photon resonances depend on the 
applied frequency v as An g = yjh 2 v 2 — Ej/8Ec, where 
we used An g = (n g hp — n peak )/2. An example of such 
an avoided level crossing is shown in Fig 2(c) where we 
have used least-squares fits to this equation to estimate 
the Coloumb and Josephson energies. 

Measurements of time-resolved, coherent oscillations of 
the charge on the SCB are made by slowly ramping the 
SCB gate charge n g as before, while applying fast (non- 
adiabatic) rectangular DC-voltage pulses of amplitude a 
in the form of a pulse train to the SCB gate. When n g 
is far away from the charge degeneracy, and a such that 
n g + a places the system at the charge degeneracy, the 
leading and trailing edges of each pulse act as successive 
Hadamard transformations. In between, the system co- 
herently oscillates between charge eigenstates for a time 
At. The trailing edge of the pulse returns the gate charge 
to n g with the system in a superposition of ground and 
excited states corresponding to the phase of the oscil- 
lation acquired during At. After the pulse, the excited 
state component decays with a relaxation time T\. The 
pulses are separated by a repetition time Tr, typically 
greater than T±. 

Since the staircase is acquired by ramping n g on a 
timescale much slower than Tr, one measures an en- 
hanced charge AQ box = Q^ x - Q° h S J x , that is time- 
averaged over the pulse train and proportional to the 
probability of finding the system in the excited state. 
Peaks in AQbox vs n g will be located at gate charges 
n peak gygjj that approximately a half-integer number of 
oscillations occur at gate charge n peak + a during the time 
At. Fig. 2(b) (bottom curve) shows a measured stair- 
case using such a pulse train. When one measures such 
staircases for varying At, then time-domain oscillations 
are evident in the At cross-sections of AQbox vs n g . This 




FIG. 3: (a) Bloch sphere showing evolution about the charge 
degeneracy for Ej /h = 6GHz, Ec/h = 9GHz and a pulse rise- 
time of 35ps. (b) Evolution for a pulse that takes the system 
past the charge degeneracy and achieves 100% charge state 
polarization (downward pointing arrow). 

is shown in Fig. 4(a), where we plot the time evolution 
near the charge degeneracy point. As Tr is increased, the 
pulse-induced charge contributes less to the time average. 
A simple calculation which ignores coherence effects re- 
maining after a time Tr leads to the dependence 

71 1 _ e-Ta/Ti 

AQ box (T R ) = 2n - i + e _ TR/Ti , (2) 

where no is the initial peak amplitude. One can estimate 
T\ at the readout gate voltage by measuring the depen- 
dence of the peak amplitude on Tr and fitting the data 
to this formula. 

Using such a pulse-train allows a measurement of the 
excited state probability even when the relaxation time 
T\ is shorter than the measurement time. Population 
averaging and mixing limit the maximum observed am- 
plitude of the AQ^x oscillations to le rather than 2e, 
and this occurs only for an ideal rectangular pulse, short 
repetition times and vanishing Ej/Ec ratio. We have 
numerically integrated the Schrodinger equation using 
pulses with a finite risetime to study the various fac- 
tors affecting both the contrast of the oscillations and 
the degree of charge state polarization. The results are 
best understood when plotted on the Bloch sphere as 
shown in Fig. 3. The dash-dotted lines show the evolu- 
tion in time of the state vector starting from an initial 
condition (upward pointing arrows) and evolving through 
a maximum polarization of the charge state (downward 
arrows). The initial state differs from pure S z due to a 
finite Ej/Ec- Fig. 3(a) shows evolution for a pulse that 
takes the system to the charge degeneracy. One finds that 
a finite pulse rise-time mimics a higher Ej /Ec ratio and 
reduces the oscillation contrast. Nonetheless, Fig. 3(b) 
shows that by pulsing past the charge degeneracy, one 
can achieve nearly 100% polarization. 

Measurements of the Tr dependence of AQbox are 
shown in Fig. 4(b), where we have used an appropri- 
ate pulse amplitude and duration to achieve nearly 100% 
initial polarization. We find T\ — 108 ± 20ns at the read- 
out point n g ~ 0.25 using a least-squares fit to Eqn.J2J). 
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The data of Fig. 4(a) shows an oscillation contrast of 
approximately 0.55e. Taking in to account the experi- 
mentally determined Ej/Ec and Ti/Tr, along with the 
measured risetime St ~ 35ps of our pulse generator (An- 
ritsu MP1763C), we find good agreement between the 
observed and expected initial oscillation contrast. Nu- 
merical solutions of Schrodinger equation show that the 
reduction of contrast due to the pulse risetime depends 
strongly on Ej. This was checked experimentally by 
varying Ej and found to be in agreement with the cal- 
culations. A fit of the data of Fig. 4(a) to an expo- 
nentially decaying sinusoid gives a oscillation frequency 
vq = 3.6GHz and decoherence time of 2.9ns. 

Using fast DC pulse trains, we were able to observe 
oscillations in time out to 10ns. The dependence of de- 
coherence time on gate charge at the top of the pulse 
is shown in Fig. 4(c) and has a sharp maximum at 
the charge degeneracy. This tells us that low frequency 
charge noise is the dominant source of dephasing in our 
qubit. As discussed and demonstrated by Vion et aZ.@], 
dephasing due to charge fluctuations is minimum at the 
charge degeneracy. The coherence time found by Vion 
et al. is much larger than that found in our measure- 
ments although their initial oscillation amplitude is much 
smaller. If we extrapolate the gate charge dependence 
of T\ (not shown) to the charge degeneracy, we expect 
T\ ~ 10ns. This implies that the decoherence time we 
measure at the charge degeneracy is dominated by the 
short relaxation time rather than pure dephasing. 

Our measured T\ was also found to be independent of 
SET bias current in the sub-gap regime. Previously, a 
considerably larger T\ was measured in a similar sample 
using a spectroscopic technique 0, and found to agree 
with the theoretical T\ ~ 1/is expected from quantum 
fluctuations of a 5011 environment. The reduced T\ ob- 
served here could be attributed to several factors. One is 
coupling to the relatively unfiltered high frequency coax- 
ial line, in addition to unwanted coupling to other metal 
traces and the microwave environment of our sample. In 
practice it is difficult to estimate the actual real part 
of the impedance on these leads at such high transition 
frequencies. Another possibility concerns quasiparticles. 
Our samples were fabricated without quasiparticle traps 
and despite the observed static 2e-periodicity, there could 
be considerable non-equilibrium quasiparticle dynamics 
occuring on these time scales. Finally, it may be that the 
use of fast DC pulses leaves the bath of charge fluctuators 
in a configuration having additional channels for relax- 
tion of the qubit. Clearly what is needed is Ti measure- 
ments using a variety of techniues — RF-rotations, spec- 
troscopy, as well as fast DC-pulse — on the same sample. 

In conclusion, we have fabricated and measured a solid 
state qubit based upon a SCB combined with a RF-SET 
readout system. Due to a relatively short Ti, continu- 
ous measurement of the SCB was employed. Fast DC 
pulses were used to coherently manipulate the qubit and 
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FIG. 4: (a) AQbox versus At. The solid line is a least-squares 
fit to an exponentially decaying sinusoid and gives a decay 
time Tb — 2.9ns and frequency vq = 3.6GHz. (b) Dependence 
of the peak height AQbox on repetition time Tr for ~ 100% 
initial polarization. The solid line is a least-squares fit to 
Eqn.(|5J giving Ti = 108 ± 20ns. (c) Decoherence time T2 
versus gate charge n g . 



time-coherent oscillations of the charge were observed. 
The initial contrast of the oscillations is relatively large 
and quantitatively understood as due to finite Ej/Eq 
combined with the finite pulse risetime. Furthermore, 
nearly 100% charge state polarization can be achieved. 
The oscillations show a maximum decoherence time at 
the charge degeneracy which indicates that charge fluc- 
tuations dominate the dephasing rate. 
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